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ABSTRACT 

Context. Massive star winds are complex radiation-hydrodynamic (sometimes magnetohydrodynamic) outflows that are propelled 
by their enormously strong luminosities. The winds are often found to be structured and variable, but can also display periodic or 
quasi-periodic behavior in a variety of wind diagnostics. 

Aims. The regular variations observed in putatively single stars, especially in UV wind lines, have often been attributed to corotating 
interaction regions (CIRs) like those seen in the solar wind. We present light curves for variable polarization from winds with CIR 
structures. 

Methods. We develop a model for a time-independent CIR based on a kinematical description. Assuming optically thin electron 
scattering, we explore the range of polarimetric light curves that result as the curvature, latitude, and number of CIRs are varied. 
Results. We find that a diverse array of variable polarizations result from an exploration of cases. The net polarization from an 
unresolved source is weighted more toward the inner radii of the wind. Given that most massive stars have relatively fast winds 
compared to their rotation speeds, CIRs tend to be conical at inner radii, transitioning to a spiral shape at a few to several stellar radii 
in the wind. 

Conclusions. Winds with a single CIR structure lead to easily identifiable polarization signatures. By contrast allowing for multiple 
CIRs, all emerging from a range of azimuth and latitude positions at the star, can yield complex polarimetric behavior. Although our 
model is based on some simplifying assumptions, it produces qualitative behavior that we expect to be robust, and this has allowed us 
to explore a wide range of CIR configurations that will prove useful for interpreting polarimetric data. 
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1. Introduction 

A corotating interaction region (CIR) arises from the interaction 
of wind flows that have different speeds: rotation of the star ul¬ 
timately leads to a collision between the different speed flows to 
produce a spiral patter n in the wind. CIRs have been well-studied 
in the solar wind (e.gjRouillard et aLbOOStfeurlaga et al.lll983l 
119971: iGosling & Pizzolll999l) . They have also been invoked to 
explain a variety of phenomena observed in massive star winds. 
Evidence of large-scale structure in the winds of massive stars 
is manifold. For example, the ubiquitousness of discrete ab¬ 
sorption c omponents (DACs) in the ultraviolet (UV) spe ctra of 
OB stars dHowarth & Prini al 1 1 989b iFullerton et al.lll997l) . some 
shown to have a re currence timescale related to the rotation rate 
dKaner et al.lll999l) . are though t to b e a co n sequ ence of the pres¬ 
ence in the wind of CIRs. ICranmer & Qwockil (1 1 996h predicted 
that such regions should develop following a perturbation, such 
as a bright or dark spot, on the stellar photosphere that propa¬ 
gates in the wind, thus generating a complex spiral-like struc¬ 
ture of low and high density and velocity regions as the star 
rotates. CIRs have also been associated with the periodic spec¬ 
troscopic, photometric, and p olarimetric variability observed in 
some Wolf- Rayet (WR) stars ( St -Louis et al.lll995l:lMorel et~aH 

119971119991:IChene & ,St-Fouisll2010h. 

Nevertheless, the origin of the perturbations generating the 
departure from spherical symmetry in massive-star winds re¬ 
mains unclear. Possible physical processes include non-radial 
pulsations and magnetic fields. The latter is particularly promis¬ 
ing. Large-scale, organized magnetic fields have been found in 


only 7 % of O stars ( Wade et all 1201 3l) . and these are thought 
to have a fossil origin. However, it has been suggested that 
smaller scale, localized ma gnetic structures could b e responsi¬ 
ble for the ons et o f DACs (Ka per & Henrichsjfl994l) . More re¬ 
cently, iCantiello & Braithwaite ( 2011 ) have presented a model 
that suggests that such structures could be widespread for stars 
that harbor a subsurface convection zone. Such zones have been 
predicted to exist in hot massive stars bv ICantiello et al.l (1 20091) 
owing to an opacity peak caused by the partial ionization of iron- 
group elements. Detecting such small-scale magnetic structures 
at the surface of massive stars is a n extremely difficult task when 
dealing with O stars dKochukhov & Sudnikl2013l) and even more 
difficult for WR stars where the photosphere is hidden by the op¬ 
tically thick wind. 


Studying the wind structures themselves can provide useful 
constraints that will shed new light on their origin. Polarimatry 
is a powerful means of studying asymmetries in stellar winds. 
Indeed, in view of their high temperatures, massive-star out¬ 
flows contain a copious number of free electrons that scatter 
light, thereby generating linear polarization. For unresolved and 
spherically symmetric envelopes, there is total cancellation of 
the polarization as integrated across the source. Consequently, 
a net polarization intrinsic to the source demon strates de viati on 
by the source geometry from the spherical. iBrown & McLean] 
( 1977 ) showed that for optically thin scattering and an axisym- 
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metric envelope, the net polarization depends on three factor^: 
the viewing inclination, an average optical depth for the scatter¬ 
ing envelope, and the shape of the envelope. 

Since electron scattering is gray, broadband continuum po¬ 
larization measurements can be used to characterize the asym¬ 
metry at any given time. If time-dependent measurements are 
available, more accurate constraints can be obtained because the 
structure will then be observed from different viewing angles. 
Spectropolarimetry can also be used to study the structure in 
the line-forming region and to obtain kinematic information. In 
the case of CIRs, this is particularly important because they are 
thought to originate in the photosphere and propagate through¬ 
out the line-forming region of the wind. 

To understand the phenomena of CIR and use their occur¬ 
rence as a diagnostic of wind and at mosp heric structure, models 
are needed to interpret observations. iMullar^ (119861) was an early 
proponent of CIRs observed in the solar wind as a viable ex¬ 
planation for certain types of variability (noted above) observed 
in massive star winds. An important breakthrough in the field 
came w ith 2D hydrodynamic simulations bv ICranmer & Owockil 
(119961) . who modeled an equator ial CIR for a ste ady-state, line- 
driven wind. Following on this, iDessartl ( 2004 ) developed 3D 
simulations for CIRs off the equator. Both of those works ex¬ 
plore the consequences for wind line-profile effe cts, su ch as 
DACs and other periodic signatures. IBrown et~akl (120041) con¬ 
sidered inverse techniques for extracting kinematic information 
about CIRs from wind emission lines. 

Pertinent to this contribution, there has also been work on the 
polarization that can arise from CIR structures. lHarries | ( 200 01) 
explored a number of effects that give rise to variable polar¬ 
ization for rotating stars, including a model for an equatorial 
CIR structure. Harries explored these effect s usi n g Mo nte Carlo 
radiative transfer simulations. In contrast, llgnace et al.l ( 20091 ) 
adopted a similar kinematic prescription for the structure of an 
equatorial CIR in an application for the variable linear polariza¬ 
tion seen in the blue supergiant star HD 92207. In this paper we 
extend these initial results in a parameter study of variable polar¬ 
ization from a wind with CIRs. Neither of the preceding papers 
allowed for multiple CIRs, and both were restricted to a CIR 
in the equatorial plane. Here results are presented that allow for 
multiple CIRs and at arbitrary latitudes. 

The structure of the paper is as follows. Section [2] presents 
the base model for the CIR structure and the resulting v ariab le 
polarization, following but expanding on llgnace etak ( 20091) . 
The results of our parameter study are presented in Section [3] 
Concluding remarks are given in Section [4] with particular con¬ 
sideration of applications for observations. Our results are pre¬ 
sented for a wind model that is set by a commonly used velocity 
structure. An Appendix provides a solution for a different veloc¬ 
ity structure as a comparison example. 


2. Polarization from a simple CIR 

In the next section, the approach o f IBrown & McLeanl (1 19771) . 
modified by ICassinelli et al.l dl987l) . for optically thin electron 
scattering (i.e., single scattering) in an axisymmetric geometry 
is reviewed. Then, an application is considered for a cone struc¬ 
ture in an otherwise spherical wind. Then we show how to use 
the cone as a basis for computing variable polarization from an 


1 The a nalysis of [B row n & McLe~anl ( 1977 1 does not include limb 

darkening. IBrown et alj d 1989h considered limb darkening, which can 
introduce chromatic effects, since it is 4-dependent. 


equatorial CIR for a rotating star, and finally the polarization for 
CIRs is explored at arbitrary latitudes. 

2.1. Polarization from optically thin electron scattering 

For optically thin scattering, the amount of scattered light from a 
circumstellar envelope is determined by integrating the volume 
emissivity over the extent of the envelope. The volume emis- 
sivity j for electron scattering depends on the electron number 
density n e , the Thomson scattering cross-section ctt, and the 
amount of incident radiation, as well as geometrical factors relat¬ 
ing to dipole scattering toward the observer. It is assumed that the 
star is the sole source of illumination for circumstellar electrons, 
meaning that the circumstellar envelope is sufficiently optically 
thin so that the amount of scattered light is small compared to 
the starlight. 

We consider a coordinate system with the Earth along the z- 
axis. The star has an axis z* that is its rotation axis. The angle 
of inclination between these axes is signified by 2 • 2* = cos ig. 
It is further assu med that y - y t , and there fore x ■ x* — cos ig as 
well. Following IBrown & McLeanl (T 977I ) for a point source of 
illumination, the scattering angle between a ray of starlight and 
the observer is^f. A scatterer is oriented about the observer’s axis 
with azimuth ijj from the x-axis. 

In the Stokes vector prescription, the polarization proper¬ 
ties o f intensity are described in terms of I, Q, U, and V (e.g., 
[Chandrasekhar ! 1960 ). The latter is for circular polarization, 
which will not be considered further. Stokes I is a measure of 
the total intensity. Stokes-Q and U describe the linear polariza¬ 
tion of the light. For example, the relative linear polarization is 


Ve 1 2 + u 2 
p = ——• 

and the polarization position angle is 


( 1 ) 


tan 2t/r p = (2) 

Again, we assume that the total flux of light will be domi¬ 
nated by direct starlight, meaning (a) that the amount of scattered 
starlight is small and (b) that the amount of absorbed starlight 
by the circumstellar envelope is miniscule. Consequently, the 
emissivity for the Stokes-I component is not a concern here. The 
emissivities for Stokes-Q and U are given by 


( ctt \ 3 . 7 

J <2 = J L v n e X - sin x cos 2 xjt, 

and 

t ctt \ 3. 2 . 

Ju = [j—zj L v”e X - sin x sin 2^. 


(3) 


(4) 


Here r is the radius for a scattering electron from the center of the 
star, and L v is the monochromatic stellar luminosity at frequency 
v. As previously noted, the angles and <// pertain to the geometry 
of the scattering. 

The star has spherical coordinates (r, i), if). The observer is 
located at cp — 0° and d - ig (see Fig. [TJ. Using g = cosi?, 
IBrown & McLeanl (119771) showed that for an axisymmetric dis¬ 
tribution of electron scatterers, integration about iA gives 

r* 'In 

I sin 2 ^ cos2i/ni(/i =-7rsin 2 /o(l - 3;U 2 ), (5) 

Jo 
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Fig. 1. Scattering geometry for Thomson scattering. Here z*, z, 
and r are unit vectors for the stellar rotation axis, observer axis, 
and local radial. The scatterer has coordinates (§, < p). The scat¬ 
tering angles are (x , <//), and i is the viewing inclination. The ap¬ 
pearance of 7t/2 ensures that a net polarization aligned with z* 
has q > 0, versus q < 0 if perpendicular. 


and 

sin 2 x sin2iA dip — 0. (6) 

These imply that only the U-polarization vanishes upon integra¬ 
tion about the scattering volume for an unresolved source, and 
only the Q-polarization survives. The total polarization p is then 
given by 

(1-3 1 ?) dr dp, (7) 

where Lq is the luminosity of polarized light, no some reference 
number density of scatterers, r - r/R t the radius in the wind as 
normalized to the stellar radius, and To an optical depth parame¬ 
ter. The factor to is given by 

t 0 = n 0 a T R t . (8) 

ICassinelli et al.l (119871) introduced a correction term to relax 
the assumption of a point source star. The correction is known 
as the finite depolarization factor, D(r). The factor assumes a 
star of uniform br i ghtnes s (i.e., no limb darkening). Although 
iBrown & McLeanl ( 119771) considered more general representa¬ 
tions for this correction, we assume a star of uniform brightness. 

The factor D accounts for the varying degree of anisotropy 
of the radiation field as perceived by a scattering electron at dif¬ 
ferent distances from the star. The factor is given by 


Lp 

~L V 


16 


t 0 sin- 1 0 


'»// 


n e (r,p) 

n Q 



D(r) = 



(9) 


This corrective modifies the polarization by being inserted as a 
“weighting” factor in the integrand of equation ([T}: 

q = To sin 2 i 0 J J j D(r) (1 - 3p 2 ) drdp. (10) 

The net relative polarization is determined by q with u = 0, ow¬ 
ing to the axisymmetry of the envelope. We note that q can be 
positive or negative, signifying the orientation of the net polar¬ 
ization against the sky. The projected symmetry axis of the en¬ 
velope is taken lie along as the north-south direction in the sky. 
Consequently, a structure like a bipolar jet would have q < 0, 
whereas one like an equatorial disk would have q > 0. The result 
shows that the polarization scales in relation to viewing inclina¬ 
tion as sin 2 i. Additionally, the polarization is linear in the optical 
depth scaling parameter for the envelope with q oc to. The ratio 
q/ro is thus a function purely from the envelope geometry. 


2.2. The case of a simple cone 

Our treatment for CIRs is in terms of segments of a cone that 
are phase-lagged from each other in a radially dependent way. It 
therefore is useful to first review the treatment of the polarization 
from a simple cone. 

We consider a cone with z * the symmetry axis of the cone: 
the cone lies directly over the stellar rotation axis with open¬ 
ing angle j3o- The star has a spherically symmetric wind, and the 
cone represents an axisymmetric sector of the wind that has ei¬ 
ther higher or lower density, thereby breaking the spherical sym¬ 
metry of the system and leading to a net polarization from elec¬ 
tron scattering of starlight as described in the previous section. In 
relation to the results presented in that section, the main consid¬ 
erations here are (a) to relate the polarization to the properties of 
the wind and (b) to use the result as a basis for our construction 
of a CIR. 

The stellar wind is taken to have mass-loss rate M and wind 
terminal speed Vco- The gas is ionized with a mean molecular 
weight per free electron p e . The number density of electrons in 
the spherical wind is given by 


«w 


M ! p e m H 
AnRlvoc 


1 


r 1 w(f) 


= no r 2 w 1 , 


( 11 ) 


where no is a parameter for the scale of the number density that 
is given by the factors for wind and star parameters at the begin¬ 
ning of the equation, and w = v(r)/v m is the normalized wind 
velocity. We use a standard wind velocity law given by 



( 12 ) 


with y > 0 the velocity law exponent, and h a dimensionless 
parameter that sets the initial wind speed at the wind base, with 
wo = (1 -b) 7 . Examples used in this paper adopt y = 1 (although 
the Appendix presents a solution with y = 2). 

The electron density in the cone is parametrized by 

n c = (1 + q)n w (r), (13) 

where q is a dimensionless parameter taking values in the inter¬ 
val q e [-1, oo). In effect, q represents an excess or decrement of 
density in the cone relative to the otherwise spherical wind. The 
cone has zero density if q = -1. If q = 0, the cone has the same 
density as the wind, in which case sphericity is preserved, and 
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one expects zero net polarization. In this form we are taking the 
run of density with radius in the cone to have the same functional 
form as in the wind. This need not be the case, and certainly a 
different radial dependence for the density in the cone as com¬ 
pared to the wind could be treated (e.g., perhaps y is different for 
the flow in the cone as compared to the wind flow). 

The results of Brown & McLean and Cassinelli et al. are gen¬ 
eral to any axisymmetric distribution of scatterers. However, a 
cone is an example of a special class of axisymmetric envelopes 
in which the density is mathematically separable in terms of ra¬ 
dial and angular functions. With p = cos 6 , we define the func¬ 
tion / and g as 

njn 0 = f(r)g(jj). (14) 

In the application here, f(r ) = r~ 2 w _1 is the same both in the 
cone and in the wind. For a cone the function g — 1 for p > 
cosflo and g = 0 forp < cos(/?o). 

Consider now the polarization from the cone. It is deter¬ 
mined by a volume integral over the entire envelope, in both the 
cone and the non-cone regions. The polarization is given by 

q = 77 t 0 sin 2 i 0 f f ” e ’ -A/> D(r) (1 - 3 p 2 ) dp dr, (15) 

16 J J n 0 

where n e is the electron density, which is either « w or n c . Since 
the cone density scales with the wind density, the integration can 
be recast as 


q = 


— r 0 sin i 0 x 
io 

r r +1 ^ 

J i J-i n 0 


D(r)(l - 3p 2 )dpdr 


r°° f +1 n 

+q I I — D{r){\-3p 2 )dpdr 

J l Ju 0 no 


(16) 


where po = cos /%. The first integral is an integration over the 
spherical wind, which of course vanishes. The second integral is 
over the cone, which scales with the factor q. 

The separability means that the integration and radius and in 
angle can be written as a product of factors, with 


q = — T o sin- i 0 F A, 
10 


where 

r = J m D{r) dr. 


and 


r 


A = g(p)(l-3p z )dp. 


(17) 

( 18 ) 

(19) 


For the problem at hand, and using a change of variable with 
£ = f~ l , these functions become 


r = 

r ~2 - lnr - w - r 

r w D(r) dr = d%, 

h Jo MO 

(20) 

and 



A = 

r +1 

1 (1 - 3 p 2 )dp = -p 0 (l - po ) 2 = A(j. 

JfM) 

(21) 



Fig. 2. Plot of the normalized cumulative polarization in terms 
of distance from the star. The figure applies when the polariza¬ 
tion calculation is seperable (see text). In the case y — 1, refer¬ 
ences lines show that half of the polarization will be determined 
within half a stellar radius above the atmosphere; 90% will be 
determined between one and 5.6 stellar radii. 


In relation to the radial integral factor, T, it is useful to consider 
the sensitivity of the polarization to structure at different radii 
in the wind. Figure [2] shows the value of F as integrated from 
£ to 1. In this form the plot displays T as the cumulative polar¬ 
ization from the stellar surface out to the radius r - £ 1 (for a 
given value of A), hence the label for the vertical axis. The plot 
is normalized to have a unit value when the integral is carried 
out to infinite distance. The reference lines show the locations 
where 50% and 90% of the polarization are obtained. In terms 
of the radial weighting, half of the polarization is set between 
1.0 and about 1.5 stellar radii; 90% is set between 1.0 and about 
5.6 stellar radii. The remaining 10% arises between 5.6 R, and 
beyond. As a result, when discussing the polarization due to a 
curved CIR, the bulk of the polarization will be determined by 
the structure of the CIR at the inner few stellar radii. 

A word about convention is useful to mention at this point. 
The observer system being used assumes the orientation of the 
cone as projected in the sky lies along a north-south direction 
(i.e., along a meridian). For a cone with q > 0, and noting that 
Aq < 0, the net polarization would be negative. This simply in¬ 
dicates that the electric vector of the net polarization is parallel 
to the east-west line, which is orthogonal to the orientation of 
the cone projected onto the sky. If q < 0, then q > 0 because the 
polarization is dominated by scatterers lying outside the cone, 
and the polarization position angle lies parallel to the cone. 

Finally, that the cone orientation is north-south corresponds 
to an orientation of the projected axis of symmetry on the sky 
as i/j - 90° as measured from west through north. Obviously, an 
arbitrary source could have an orientation on the sky, and the net 
polarization would not change. It is straightforward to represent 
this polarization and position angle as follows. Let po - \q(i// - 
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90°)|. Then for arbitrary orientation as perceived by an observer, 
the normalized Stokes q and u parameters become 


q = po cos(2i//), and ( 22 ) 

u = po sin(2t/r). (23) 

In developing a simplified model for the polarization from CIRs, 
we use these constructions in what follows. 

2.3. CIRs in the equatorial plane 

A CIR amounts to a spiral-shaped zone of flow that threads a 
stellar wind. Its structure can inde ed be complicated, as seen 
in t he hydrodynam ic simulations of ICranmer & Owockil (1 19961) 
and iDessarn (12004 ). Our goal is to construct an approximation 
for the structure of a CIR to rapidly explore model parameter 
space and trends in the variable polarization with rotation phase. 
To this end, we approximate a cross-section of a CIR with a 
sphere as a “spherical cap”. This cap has circular cross-section 
and is identical to the intersection of a sphere of radius r with a 
cone. Our construction for a CIR is simply these cap segments 
for a cone, but where each successive segments in radius are az- 
imuthally phase-lagged. In fact we require an equation of motion 
for the centers of the segments. 

Initially, we consider a CIR that is spiral in shape where the 
center of the spiral lies in the equatorial plane. The calculation 
for the polarization can take an arbitrary function of density into 
account with radius inside the cone and even a varying solid an¬ 
gle of the segments with radius. As long as the CIR segments 
are individually axisymmetric and scattering is optically thin, 
the approach described in preceding sections can be employed. 
However, to reduce the number of free parameters, we assume 
that all segments of the CIR have a constant solid angle, and we 
adopt the same parameterization for the CIR density as for the 
cone of the preceding section: the scaling with q such that 77 is a 
constant throughout the CIR, and the y velocity exponent is also 
the same in both the cone and the wind. 

The question now is what to use for the equation of motion 
of the segment centers. What is the geometric form for the spi¬ 
ral? A physically m otivate d approach could be based on wind- 
compression theory dBiorkman & Cassinellilfl993t llgnace et al.l 
1 19961 ). That model is for a rotating wind in which streamline flow 
is solved via kinematic relations, namely a velocity law and con¬ 
servation of angular momentum. This model for rotating winds 
leads to density enhancements toward equatorial latitudes at 
the expense of polar latitudes. However, wind-compression the¬ 
ory assumes a central force for the wind driving. lOwocki et al.l 
( 1996 ) have shown that this does not hold for line-driven winds 
for massive stars. Indeed, Owocki et al. find that including non- 
radial forces leads to density-enhanced polar flows. 

For our purposes we adopt purely radial streamlines for wind 
flow in the observer’s frame. This means that a parcel of fluid 
injected into the wind from the star travels along a radial line, 
but since the star rotates, subsequent parcels of gas also follow 
a radial trajectory, but are successively displaced in azimuth. 
Consequently, the CIR amounts to a spiral pattern of flow that 
spins as a whole with the rotation period of the star. It is basi¬ 
cally the same problem as for perceiving a spiral pattern of water 
flow from a rotating hose. The shape is found from the stream¬ 
line flow for a CIR in the rotating frame of the star and recog¬ 
nizing that the shape is equivalent to the spiral pattern as se en 
by an observer. (A similar approach was used bv llgnace et al.l 


( 1998 3 to model the magnetic field topology in the context of 
wind-compression theory.) 

In the rotating frame of the star, we imagine a generalized 
“spot” of opening angle /?o and solid angle ilo at the equator of a 
star that rotates with period P and angular velocity to. This spot 
is the source of a density perturbation on the wind (i.e., the cause 
of q 4 0). In the rotating frame, the spot is fixed. If we assume 
that the flow from the spot only has a radial component, then in 
the rotating frame, a fluid element emerging from the spot will 
progressively displace the spot azimuthally in proportion to tut, 
where t is the interval of time after the fluid parcel was injected 
into the flow. 

Generally, the fluid element should of course emerge from 
the spot with an azimuthal speed of tuR „. However, with the com¬ 
plexity of non-radial force considerations, we simply ignore the 
details of the azimuthal component of speed v^ that a fluid el¬ 
ement would have. Certainly, a parametrization for \p could be 
included. In fact, for v, f the radius-dependent azimuthal velocity 
in the frame of the star, the equation of motion for the center of 
the fluid element in the rotating frame of the star would be 

dip' 1 r vy(r)-a>r l 

dr r v r 


where <p' is the azimuth about the star in the rotating frame as 
measured from the spot (i.e., ip' - 0 ° is the location of the spot). 
Our model for the spiral CIR assumes v v (r) = 0. 

The solution to the preceding equation then becomes 


r 

<p 



1 


e w© 




(25) 


where ro — Voo/to is a convenient parameterization representing 
the “winding radius” of the CIR spiral. Fast radial flow means 
that the spiral is mostly a linear cone at radii close to the star be¬ 
cause the rotation is relatively slow, and so ro/R* is large, yield¬ 
ing ip’ is roughly constant for r < ro. In contrast, fast rotation 
relative to the radial flow speed indicates a rapid winding up of 
the spiral, which is represented by a relatively low value of ro. 

In the observer’s frame, the equation for the center of the 
spiral is 


ip - tpo + tp'(r) + tut. 


(26) 


where ipo allows for the azimuthal location of the spot to be de¬ 
fined in terms of an arbitrary reference. 

For -y=l, the integral equation in £ is analytic. The solution 
is 


ip'(r = £ *) = - — 

ro 


w 


+ b In 



(27) 


Analytic solutions can also be found for other integer values of 
7 (see the Appendix). For our case of 7 = 1, the influence of the 
winding radius, ro/R,, for the shape of the CIR spiral is shown 
in Figure [3] Two reference circles are shown for radii of 2 R, 
and 4 R,. The four curves are for the solution to tp' for different 
values of ro/R * = 2,4,6, and 8 , as labeled. Lower values of ro 
indicate relatively high equatorial rotation speeds as compared 
to the wind terminal speed. Curves with lower ro have higher 
values of ip' at a given radius. 

To evaluate the polarization, a key result is to recognize the 
separability of the spatial integrations for our adopted model. 
Our assumption is that the cross-section for the CIR is the same 
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Fig. 3. Examples of the spiral shape for the CIRs as a function of 
the winding radius, r^. The figure displays the centroid of four 
different CIRs with distance from the star. The star is the shaded 
magenta circle. The CIRs are in the equatorial plane, shown here 
as seen from above. The two larger dotted circles are reference 
circles at 2R t and 4 R t . The different spirals are distinguished by 
different winding radii as labeled. Lower values of r (l place the 
transition from a cone-shaped CIR to a spiral one closer to the 
star; higher values of ro indicate that the CIR has small curvature 
until achieving large distances from the star. 


as for a cone; therefore, Ao for the cone holds for the CIR. The 
polarization for a CIR reduces to evaluating T. However, the so¬ 
lution for T, which is an integration in radius, must take account 
of the radius-dependent phase-lagging of segments as given by 
the solution of equation (l27l >. 

The phase-lagging must be related to the observer system 
of coordinates. The result from Brown & McLean reveals that 
the polarization from a given conical segment scales with sin 2 i. 
For a CIR, each segment is essentially at a different inclination 
and position angle. In effect, each segment has i = i(ip(r, t )) and 
i fi = t)). Using spherical trigonometry, one can show that 


sin 2 /cos2i/r = sin 2 ip - cos 2 i'o cos 2 ip, and (28) 

sin 2 i sin 2if/ — - cos t'o sin 2i p. (29) 

Then the polarization for a CIR at a given a time in terms of the 
normalized Stokes parameters q and u is given by 


<7(0 = 

qT Q A 0 T q , and 

(30) 

u(t) = 

77 T 0 A 0 T U , 
lo 

(31) 

where 



r q - _ 

f 1 £>(£) , 

I - sin" i cos 2i // d£, and 

Jo 

(32) 


T u = f —— sin 2 / sin 2 if/dt;, (33) 

Jo w(f) 

with 

sin 2 /cos2i/r = sin 2 ip(t, g) - cos 2 io cos 2 ip(t, f), and (34) 
sin 2 /sin2i/r = -cos/o sin2ip(t, £). (35) 

The total polarization and the polarization position angle become 

P = JffP T o A <) >/r 2 +T 2 , (36) 

and 

tan 2 i/r = - = (37) 

q r q 

where the last equality is specifically true for the many assump¬ 
tions that we have adopted, such as every cross-sectional seg¬ 
ment of the CIR having the same solid angle. In the optically thin 
limit, the polarization position angle is independent of the wind 
optical depth, the value of q, and even the CIR cross-section as 
encapsulated in Ao. The only relevant factor is the shape of the 
spiral as a function of radius as weighted by the geometrical fac¬ 
tors associated with Thomson scattering. 

It is straightforward to include multiple CIR structures by 
computing the q and u values for each one separately. This holds 
as long as each one is optically thin. For N CIRs at the equator, 
each one having its own value of //,, i^o.i, and fi, for the i' h CIR, 
the Stokes q and u parameters are linear combinations with 


*7tot 


^tot 



(38) 

(39) 


i=l 


However, it is useful to formulate the total mass loss of the flow. 
The angle-averaged electron optical depth of the entire envelope 
becomes 


T e = T„ 




i=l 


(40) 


If one associates a certain M with the wind component, then a 
total mass-loss rate can be determined from the preceding equa¬ 
tion after inclusion of CIRs. 

It is worth mentioning that CIRs arise from redistribution of 
wind flow, as contrasted to a model of sectors with enhanced or 
diminished mass loss relative to an average spherical wind. In 
redistribution one envisions CIRs as perturbations of the wind 
density, and the summation term in equation (l40l) must vanish. 
In this framework a “CIR” would actually consist of two CIR 
structures, one with positive q and the other with negative q so 
as to conserve the mass-loss rate. In other words redistribution of 
density implies that f e = r w , and our model would then suggest 
pairs of CIRs with high {q\f) and low (q\) densities, such that 

/7 h C2h + q\£l\ = 0 (41) 

for every pair. Model results to be presented in section 3 adopt a 
framework of high and low density streams instead of a consid¬ 
eration of redistribution. The point is that redistributive approach 
could be modeled in a crude way with pairs of CIRs as just de¬ 
scribed. 
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2.4. CIRs from arbitrary latitudes 

Off the equator, the processes for formulating the CIR geometry 
and the polarization is largely the same. At the equator the CIR 
remains centered in the equatorial plane. At other latitudes the 
center of a CIR moves along a conical surface under our scheme 
(since the flow is radial). The equation of motion for ip(t, r, Sq) 
for a stream originating at latitude So is almost unchanged from 
the equatorial case. With nr = r sin So, the differential equation 
for the centers of CIR segments is 

d -f = --- (42) 

dm v r 

Given that for our approximation, 6 = Oq is fixed at all radii for 
a CIR originating at So, the solution for p becomes 

R CcU: 

<p(t, r. So) = (Jjt + (fio - 1 Sin 8 0 -r—. (43) 

ro Jt £ w 

The main difference is the appearance of the factor sin Sq. The 
interpretation is that there is now an “effective winding radius” 
with r e ff = ro/ sin Sq, and this becomes larger for latitudes closer 
to the pole. At the pole, r e ff —> oo, and a CIR emerging from the 
pole would be a conical stream. 

A major difference from the equatorial case comes in the 
form of the more complex spherical trigonometric relations that 
transform between the star coordinate system and that of the ob¬ 
server. The expressions needed for T q and T u become 

sin 2 i cos 2>p = sin 2 So sin 2 p - (cos So sin io 

- cos io sin So cos ip) 2 , and (44) 

sin 2 i sin 2 ip = sin io sin 28q sin ip 

- cos i'o sin 2 So sin 2tp. (45) 

With these relations in hand, polarization light curves can be 
computed for a CIR emerging from an arbitrary location on 
the star. Indeed, using the approach of the preceding section, a 
model with multiple CIRs can also be calculated. 

The next section details an exploration of parameters as a 
sampling of the kinds of variable polarizations that can result for 
a number of select scenarios. 

3. Results 

We have conducted a broad parameter study for the variable po¬ 
larization resulting from our model of CIRs. Before discussing 
these results, we wish to emphasize two points. First, these 
models represent the intrinsic source polarization. Real mea¬ 
surements are contaminated by an additional component aris¬ 
ing from interstellar polarization. An unresolved star that is 
spherically symmetric will generally have a measurable polar¬ 
ization of the polarizing influence of the interstellar gas and dust 
through which the starlight must travel to reach Earth. This influ¬ 
ence depends on wavelength. The effec t is frequently modeled 
or removed using a “Serkowski Law” (ISerkowski et al.lll975l) . 
Suppose the star has an intrinsic normalized polarization given 
by [qJA), uJA)] at the Earth. The interstellar contribution can be 
modeled as adding an additional [qisM('l), m/sm( 4)]. The mea¬ 
sured polarization will then become ( q » + qjsM, u* + ujsm )■ The 
ISM contribution essentially produces a wavelength-dependent 
translation of the intrinsic stellar polarization in the q - u dia¬ 
gram. 



Phase q/ T() (%) 


Fig. 4. Periodic variations of polarization with stellar rotational 
phase for CIRs with different winding radius values. Left panels 
are for polarization light curves with phase; right panels are for 
variations in the q-u plane. Top panels are for a viewing inclina¬ 
tion of i = 30° and lower for i — 60°. Different curves are for 
different winding radii, with values ro/R* — 1 (black), 3 (red), 9 
(green), and 27 (blue). 


Second, and most important for our purposes, the interstellar 
contribution is constant. Any variable polarization is therefore 
necessarily intrinsic to the stellar source. One cannot generally 
say that the interstellar polarization is given by the observed min¬ 
imum q and u values. In other words, even if the star did not 
display variable polarization, it may still have a non-zero intrin¬ 
sic polarization. Nonetheless, variable polarization is taken as 
evidence for intrinsic polarization. If the amplitude of that polar¬ 
ization is larger than the minimum value, then one may certainly 
conclude that the stellar polarization is large compared to the 
interstellar component. 

With these points in mind, what follows presents idealized 
model results for q jt) and ujt) that ignore the ISM component, 
as well as considerations of sources of random or systematic 
(e.g., instrumental) errors present in the data. Our assumption 
is that our models would be applied to processed data after in¬ 
terstellar and systematic effects have been corrected. 

The parameter space for CIR simulations is large owing to 
our rather general model. Our model for the structure of CIRs is 
admittedly parameteric. We have strived to minimize the num¬ 
ber of free parameters in our approach. Nonetheless, simulation 
variables include opening angle /?, viewing inclination i, latitude 
of the CIR §, azimuthal location tp, the winding radius ro/R t , 
and the density enhancement factor //. In addition, simulations 
allow for an arbitrary number of CIR structures N. 

As a sampling of this enormous range of possible models, 
we focus on four basic sets: the effect of changing the winding 
radius for one equatorial CIR, the effect of multiple CIRs all 
emerging from the equator of the star, the effect of CIRs from 
different latitudes as seen from different viewing inclinations, 
and the effect of different numbers of randomly distributed CIRs 
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Fig. 5. Periodic variable polarization for two equatorial CIRs. 
The format is the same as in Fig. [4] with light curves to the left, 
q-u plane variations to the right, a viewing perspective of i = 30° 
at the top, and i - 60° along the bottom. The two CIRs are 
shifted from each other by an angle A tp of 0° (black), 30° (red), 
60° (green), 90° (dark blue), 120° (light blue), 150° (magenta), 
and 180° (orange). 


as seen from a typical viewing perspective. Except for the first 
set of models, the winding radius will be fixed at nj R r = 5. 
This amounts to an equatorial rotation speed that is 20% of the 
wind terminal speed. In addition, the opening angle [5 = 15° is 
fixed. Thus at a given radius, the CIR covers a solid angle of just 
under 2% of the “sky” as seen from the star. Since our model 
is for optically thin scattering, the polarization parameters q, u, 
and p all scale linearly with the optical depth to, so results will 
be presented as normalized to to as percentage values. 

3.1. Changing the winding radius 

Figure [4] shows model results for a single equatorial CIR as the 
wind radius ro/R» is varied. Left panels are polarization light 
curves; right are q-u planes. The top panels are for a viewing 
inclination of i = 30°, whereas the bottom panels are for i — 90°. 
The curves are for ro/R t - 1,3,9, andll. Higher values for the 
winding radius imply that CIRs start to curve at larger radii. 

Not surprisingly, small winding radii tend to have low rela¬ 
tive polarization (i.e., p/tq) amplitudes. This is because material 
is more spread out in azimuth about the star, and so a greater 
range of scattering angles are sampled, leading to polarimetric 
cancellation. Higher values of ro/R* lead to higher relative po¬ 
larization amplitudes. In all cases, two peaks separated by half a 
rotational phase are seen in the polarization light curve because 
there is a single CIR structure. In the q-u plane, these lead to a 
double loop. The two troughs of the light curves (and the two 
loops in the q-u plane) are not identical because of stellar occul- 
tation in conjunction with the CIR being asymmetric. 



Phase q/ T() (%) 


Fig. 6. Periodic variations in polarization are shown here for 
CIRs that emerge from different latitudes of the star. The layout 
of the figure is the same as in Fig. [5] Here the colors are for CIRs 
at latitudes of i? = 20° (black), 40° (red), 60° (green), and 80° 
(blue). 


3.2. Equatorial CIRs 

We consider a model with two equatorial CIRs. The first CIR has 
q — 1 and emerges from an azimuth of - 0°. A second also 
has q - 1 but emerges from (p — 0° to 180° in 30° increments. 
(Models from 180° to 360° are degenerate with our set, within 
a phase shift.) The two are co-added with q - q\ + qi and u = 

Ml + M2. 

Figure 0 displays polarimetric light curves for this set of 
models. The display is the same as in Figure[4] The different col¬ 
ors for the different azimuthal offsets between the pair of CIRs 
is indicated in the caption. All of the light curves are double 
peaked, but (a) the shape of the peaks, (b) the interval between 
the peaks in phase, and (c) the amplitude of the peaks depend on 
the relative positioning of the two CIRs in azimuth. Of particu¬ 
lar interest is that for the parameters shown, the two peaks are 
closest when the CIRs emerge at 90° from one another around 
the equator. 

3.3. CIRs at different latitudes 

In this set of models, a single CIR structure threads the wind, but 
with the CIR emerging from different latitudes and viewed from 
different inclinations. Figure [6] shows results in a format similar 
to Figure [4] Like that figure, the upper panels are for i = 30° 
and the lower for i = 60°: left is for polarization light curves and 
right for the q-u plane. The different curves are for CIRs from 
different latitudes, with d = 20°, 40°, 60°, and 80°. 

Changing the latitude has distinctive effects on the variable 
polarization as compared to the previous two sections that con¬ 
sidered only equatorial CIRs. Changing the latitude, all else be¬ 
ing equal, has consequences for whether the light curve is single 
or double peaked. In the q-u plane, the effects are either a single 
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or a double loop. The reason is that a CIR that is more closely 
located near the pole tends to be more nearly stationary in pro¬ 
jection on the sky, despite the rotation of the star. The limiting 
case is a CIR at the pole, which is just a conical flow that main¬ 
tains a fixed projected position angle. 


3.4. Random distributions of CIRs 

The case that may perhaps most closely represent real stars is 
one with multiple CIR structures emerging from different az¬ 
imuth and latitude locations. In this example seven models have 
been calculated with individual CIRs. A random number genera¬ 
tor was used to place the individual CIRs around the star. Values 
of tj were randomly assigned in the range of -1 to +4. Table Q] 
summarizes the assigned parameters for the seven cases. 

Figures[7]and[8]show polarization light curves and q-u planes 
using the following construction. The viewing perspective is 
fixed at i = 60°. All the CIRs have ro/R* = 5. Then the results 
are superimposed and normalized according to 


*7 net 



(46) 



and 


_ 1 

^net — — 

n 


I 

/= 1 


U[ 


(47) 


where n runs from 1 to 7. Figures [7] and [8] thus show results 
for one CIR, two CIRs, three CIRs, and so on, up to all seven 
CIRs together. This is an incremental approach for seeing what 
happens as additional randomly placed CIR structures are added 
to the model. 

The normalization ensures that the curves are comparing the 
same things. For example, if we were to co-add seven models of 
exactly the same CIR, the resulting polarization would simply 
be the same curve as one CIR, only with a polarization that is 
seven times greater. Normalization ensures that as more curves 
are added, the curves are being compared on the same basis. 

The result is a diverse set of possible outcomes. The black 
curve shows the case of just one CIR. Adding additional CIRs 
moves through the curves in color in the order of red, green, 
dark blue, light blue, magenta, and finally orange. A variety of 
shapes with more or fewer peaks result. Add to this the fact that 
for this example, all CIRs have the same opening angle and that 
only one viewing inclination case is being shown, it seems clear 
that a wide variety of variable polarizations could result. 


4. Conclusion 

In this paper we have described a parametric representation of 
a corotating interaction region (CIR) that threads an otherwise 


Table 1. Model parameters for randomly distributed CIRs 


Case 

rj 

V (degs) 

cos & 

a 

-0.64 

220 

+0.18 

b 

+ 1.83 

250 

-0.89 

c 

+0.20 

219 

-0.22 

d 

+0.48 

290 

+0.98 

e 

+0.63 

14 

+0.22 

f 

-1.94 

326 

-0.63 

g 

+0.52 

154 

-0.74 


Fig. 7. Periodic polarization light curves for a randomized distri¬ 
bution of CIRs across the star. Each curve represents the addition 
of one more randomly placed CIR. The curves correspond to a 
number of CIRs of 1 (black), 2 (red), 3 (green), 4 (dark blue), 
5 (light blue), 6 (magenta), and 7 (orange). Table [I] lists model 
parameters of the individual CIRs. 


spherical wind. There is no denying that we took some liber¬ 
ties in presenting a simplified structure for a CIR: circular cross- 
section, constant solid angle, uniform density within the CIR, 
and a simple prescription for its spiral curvature. However, the 
advantages of the approach allow, for the first time, a broad 
exploration of consequences for variable polarization signals. 
The goals were to adopt a flexible, semi-analytic description 
of a CIR structure, basically motivated by hydrodynamic sim¬ 
ulations, in which to explore the range of polarimetric be¬ 
havior easily and rapidly. These results may then be used in 
conjunction with other methods, suc h as U V line variability 
(for a r ecent example, see iPrinia et alJl2012h . optical variabil- 
ity (e.g., Chene_&SfHouiji20j 0|), and perhaps X-ray variability 
(llgnace et al.ll20l 3^ Massa et al.ll2014h . 


An important result of our study is allowance for CIRs from 
latitudes other than the equator. The majority of previous works 
have focused on equatorial CIRs (one exception being iDessartl 
120041) . The motivation is often one of simplicity. Frequently, 
CIRs are considered in relation to variable blueshifted abso- 
prtion for UV line data involving the intersection of the CIR 
with the absorption column between the observer and the star. 
Consideration of an equatorial CIR makes the modeling eas¬ 
ier. In contrast, polarization is sensitive to the fully three- 
dimensional structure of the asymmetric distribution of density 
about the star. (The trade-off in relation to line-profile studies is 
the loss of velocity shift information.) A major consequence of 
having CIRs from different latitudes rather than just the equator 
is a far more diverse range of polarimetric light curves and of 
behavior in the q-u plane. The potential-added complexity can 
be a strength when used in conjunction with other diagnostics, 
such as line variability. 
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q/ T o ( % ) r / R* 


Fig. 8. Polarization model results from Fig. [7] displayed in the 
q-u plane. The colors correspond to the identifications indicated 
in Fig. [7] 


This contribution focused on the theoretical aspects of CIRs 
and variable polarization. There are several data sets to which 
these results could be applied. In a follow-up paper, our model 
will be employed to interpret variab le pol arizati on for the Wolf- 
Ray et (WR) star EZ CMa (iDrissen et al.lll989h . WR stars are 
known to form pseudo-photospheres, where optical depth unity 
occurs in the wind itself fe.g JHillieri 19871) . Although our model 
is for optically thin scattering, useful results can still be obtained 
through adoption of the “last scattering approximation” (e.g., 
lStenfloil982b . The merit for such an approach is empirically sug¬ 
gested bv ISchulte-Ladbeck et al.l (119921) for EZ CMa. However, 
the approach will require modifications to the prescription for 
CIR curvature, which will be explained in the follow-up paper. 
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Appendix A: CIRs for the Case y - 2 


Our results have focused exclusively on the case of y = 1 for 
the wind velocity law. Of course, other values of y may be con¬ 
sidered. In general, the solution for the spiral shape of the spiral 
will not be analytic but must be evaluated numerically. However, 
analytic solutions for < p{^) can be found for integer values of 
y. Expected values of y f or early-type star winds range from 
y = 0.5 from the original ICastor et ah (|l975 l) pape r for line- 
driven winds to y « 0.8 from iPauldrach et al.1 d 1986b that aug¬ 
mented the initial results of Castor et al. to the y = 1 tha t is typ- 
ically used for the inner wind of WR stars (e.g.. ISchmutz et aD 
119891) and then to y ~ 3, which has be en suggested for some 
supergiant winds (e.g., IPrinia et ahl 19951) . 

As an example, the case of y = 2 is provided here as a com¬ 
parison case to y = 1. In general, higher values of y tend to in¬ 
crease the radial width over which the bulk of wind acceleration 


Fig. A.l. Comparison of the azimuth location of the CIR cen¬ 
troid with radius from the star. Here blue is for y = 1, and red is 
for y = 2. This example uses a wind radius of r (l = 1 OR.. The 
value of 4>' is generally between 50% to 100% larger for y = 2 
than for y = 1 indicating that a spiral CIR is considerably more 
curved at a given radius for the higher y case. 



Rotational Phase 

Fig. A.2. Comparison of polarization light curves for the y = 1 
and y = 2 cases shown in the previous figure. The polarization 
is plotted as normalized to the wind optical depth along a radial. 
The case of y = 1 produces a relatively stronger polarization 
signal (per unit optical depth) than in the case of y = 2. The 
polarization extrema are slightly shifted in phase between the 
two cases. 
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takes place. For example, with w — (1 - b/r) r , the radius in the 
wind where the flow achieves half its terminal speed, w = 0.5, 
occurs at f\ji — b/( 1 - 0.5 ^ y ), which increases as y becomes 
larger. As a case in point, f \/2 — 2 b for y — 1, but for y — 2, the 
value of f ] i 2 increases to 3.4 b. The net effect of this is that for 
our construction of the CIR spiral shape, at a given value of a>, a 
CIR tends to wind up at smaller radius with increasing y, so the 
effective winding radius decreases with larger y. 

We need to evaluate <p' for y — 2. We start with Eq. (??) for 
a CIR in the equatorial plane, reproduced here: 



The wind velocity is w — (1 - bt;) 2 . Inserting w into the above 
expression gives 


f' 1 

* r o J( g 2 (l-bg) 2 

The solution for this expression is 


(A.2) 


/ 

<P 


/?, 17 2b - 1 \ 
ro [\ 1 - b I 


1 -2bf \ / l-b£ 

+ 2 b In f 


f(l -H) 


fd -b) 


(A.3) 


Figure lATl shows a comparison between these functions. There a 
value of ro//?* = 10 is used. At most radii the value of <p' is about 
50% to 100% higher for y — 2 than for y — 1. As compared to 
7 = 1 , the effective value of ro is about half as large when y — 2 . 

For the polarization one should also note that for a given 
optical depth of the envelope, larger y essentially implies a rel¬ 
atively higher density of scatterers at small radii as compared to 
winds with lower y values. Figure IA.2I compares the polarization 
light curves for the solutions displayed in Figure lATl equatorial 
CIRs, ro//?* = 10, wo = 0.03, same opening angle, with 77 = 1, 
and i(j = 60° for both cases. The polarization is plotted as the ra¬ 
tio p/t, where r is the optical depth of the wind. As can be seen, 
for the selected parameters that are relevant to fast winds, the 
polarization per unit optical depth is similar for the two cases, 
but there is a slight phase shift owing to the different degrees of 
CIR winding. This result does not assume the same base density 
no for the two cases, but rather the same wind optical depth. (For 
the same no, the optical of a 7 = 2 wind is 1.6 times greater than 
for 7 = 1 .) 
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